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Conduction in bulk poly crystalline high-T c superconductors with relatively high critical currents has 
been shown to be percolative. This phenomenon is due to weak links at grain boundaries. These 
weak links are the major limiting factor for technological applications which require high current 
densities. We formulate a model of these materials which can be reduced to a nonlinear resistor 
network. The model is solved by analytical approximations and a new numerical technique. The 
numerical technique is variational, which makes it capable of solving a wide variety of nonlinear 
problems. The results show that the presence of a distribution of critical currents in the sample 
does not erase all information about the dissipative electrical properties of individual boundaries. 
This means that an unambiguous connection can be made between the I — V characteristics at the 
microscopic level and the macroscopic electrical properties. 



PACS Nos. 74.25.Fy, 74.62.Bf, 74.60.Jg, 74.72.-h 



I. INTRODUCTION 

The discovery of high-T c superconductivity raised 
hopes of an important role for these materials in applica- 
tions requiring high current densities, such as high- field 
magnets and transmission lines. The most important 
factor limiting current densities in present-day low-T c 
materials is flux motion, and the solution is to maxi- 
mize pinning. While this is important in the high-T c 
materials as well, an even more vexing problem is that 
of weak links. Practical materials are polycrystalline: 
grain boundaries and other extended defects are unavoid- 
able. These boundaries, usually not important in low- 
Tc materials, limit the flow of supercurrent in high-T c 
systems. While bulk samples may be superconducting 
at low current densities, their critical currents are low. 
Since a transport current must pass through many of 
these boundaries, they pose a crucial problem in conduc- 
tor development. 

Experimentally, early work on grain boundaries 
showed that their critical currents I c could vary over or- 
ders of magnitude Gj. Under these conditions, it might 
be expected that supercurrent flow in polycrystalline 
samples, at least near the critical current, would be per- 
colative in nature. Recent work using magneto-optic 
methods has confirmed this. Particularly striking is work 
on BSCCO/Ag tapes Q|. As a result, one must under- 
stand the behavior of current flow when a distribution of 
boundaries is present, with possibly a wide spread in I c . 

Onthe theoretical side, percolation theory has dealt 
mainly with linear circuit elements and focused on the 
critical behavior (i.e., power laws) near the percolation 
threshold ||. Some work has been done with nonlinear 
elements. An example is the case of a lattice of resistive 



elements with arbitrary nonlinearity [B, || where indi- 
vidual resistive elements are removed at random. The I- 
V characteristic can be shown via renormalization group 
theory to follow a power law behavior near the percola- 
tion threshold. These renormalization group calculations 
do not tell us anything about the value of the critical cur- 
rent, or the shape of the I — V characteristic away from 
the critical region. Some work away from threshold has 
been done for metal- insulator problems ||. This prob- 
lem is in a certain sense dual to the superconductor-metal 
problem studied in this paper, but the exact relation is 
not clear. The case of a linear medium containing a small 
admixture of nonlinear elements has also been studied. 
This model has applications in the study of composites 
formed of nonlinear impurities embedded in a linear host 
1 @1 ) H ■ Theoretical predictions of I — V characteris- 
tics are made using the Clausius-Mossotti approximation 
0, or the effective medium approximation p|. These are 
limited to the case of weak nonlinearity or a low density 
of nonlinear impurities. 

The critical phenomena approach and the weak non- 
linearity case are not closely relevant to our goal here, 
which is to understand the entire I —V characteristic of 
a strongly nonlinear system. More closely related is the 
work of Leath and Tang. These authors first considered a 
Ginzburg-Landau model of conduction in polycrystalline 
superconductors, and later a simpler nonlinear resistor 
network |nj. The latter have generally been the more 
popular since then, since they provide a reasonable and 
relatively simple representation of the physical proper- 
ties of a superconductor with weak links. Hinrichsen et 
al. considered a model which is a special case of the ones 
we shall treat [jy], and it will be discussed further below. 
Our goal is different from both these papers, however. It 
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lies in determining the connection of the "microscopic" 
distribution of link strengths to the macroscopic I — V 
characteristic. That is, we are interested not in a partic- 
ular underlying model, but rather in determining what 
that model is, given global information. Hence we look 
for solution methods which will apply to whole classes of 
models, and which therefore amount to generalizations 
of those which have been discovered to date. 



II. THE MODEL 

The current flows through a random array of micro- 
crystallites separated by grain boundaries which act as 
weak links. In this paper, we shall restrict our attention 
to the two-dimensional case - a granular superconducting 
film. There is nothing which prevents the generalization 
of any of the arguments to higher dimensions. The com- 
putational method we use is independent of dimension- 
ality. 

The first step is to construct an equivalent resistor net- 
work. For the purposes of simplicity we shall assume 
that each grain is perfectly conducting. This allows one 
to separate, at least temporarily, the two current-limiting 
factors (flux pinning and weak links) mentioned above. It 
means that our results are limited to cases in which the 
total magnetic field (applied plus self) is small enough 
that there is no flux penetrating into the grains. Flux 
only penetrates the boundaries. 




FIG. 1. Reduction of a poly crystalline superconducting 
network to an equivalent resistor network. The irregularly 
shaped regions are superconducting grains, taken as equipo- 
tentials. All potential drops take place across grain bound- 
aries, shown as (nonlinear) resistors. 



Once the grains are taken as perfectly conducting, then 
they are equipotential areas. They are the nodes of the 
network. Each pair of grains is separated by a boundary. 
Each boundary is replaced by a (nonlinear) resistor. An 
example of this construction is shown in Fig. [lj. In graph- 
theoretical terms, the graph of boundaries is replaced 
by its dual graph. A theorem from topology states 
that this relation is one-to-one. We shall suppose that 
the geometrical randomness is less important than the 
randomness in the strength of the resistors. Accordingly, 
we shall consider only square lattice networks. 

Although we will develop techniques applicable to re- 
sistive elements of arbitrary nonlinearity, our main inter- 
est in this paper is the flow of current through a gran- 
ular superconductor. Accordingly, we focus in this pa- 
per on two models of the grain boundaries applicable to 
the cases of Josephson junctions in low and high applied 
fields. Why actual boundaries should behave in these 
ways is discussed by Likharev jl6| . 




Current 

FIG. 2. I — V characteristics for (A) Josephson junction 
(JJ) and (B) flux flow (FF) models. The former is considered 
to be appropriate for low fields, and the latter for high fields. 
In the numerical calculations, we give the JJ characteristic a 
finite slope just above the critical current for numerical sta- 
bility. 

Each boundary is specified by its I — V characteris- 
tic. V is some nonlinear function of I for each resistor, 
but each resistor may be different. In the general case, 
the resistors are drawn from a probability distribution 
P[V(/)]. We will concentrate on binary distributions, 
since the results easiest to understand. These binary dis- 
tributions consist of a mixture of weak links and strong 
links, each weak link being identical, and each strong link 
being identical. The concentration of weak links will be 
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denoted by p and the concentration of strong links by 
q = 1 -p. 

We further divide our model into two classes cor- 
responding to weak-field Josephson junction (JJ) and 
strong-field flux flow (FF) boundaries. For the JJ case, 
the weak links, with concentration p, have a critical cur- 
rent Id- If I < Id in such a resistor, then V = 0. If 
/ > 7 c i, then V = IRi, where Ri is a constant. Simi- 
larly, the strong links, with concentration 1— p, satisfy: if 
/ < I c2 , then V = 0; if I > I c2 , then V = IR 2 - I C 2 > Id- 
Nate that V is a discontinuous function of I for both 
types of links. Let us call this the binary JJ model. For 
the FF case, the weak links, with concentration p, have 
a critical current I c \. For I < 7 c i, V = 0. For I > 7 c i, 

V = (7 — Id)Ri> where R\ is a constant. The strong 
links, with a concentration 1 — p, satisfy: if / < I C 2, then 

V = 0; if / > I c2 , then V = (7 - 7 c2 )i? 2 . 7 c2 > 7 cl . 
Note that V is a continuous function of / for both types 
of links. This is the binary FF model. See Fig. || for an 
illustration of the two cases. 

The two models represent very different microscopic 
pictures. The JJ model represents a low-current, zero 
applied field situation. The junction is taken as the sim- 
plest sort of Josephson junction. It is superconducting 
for small currents. It is normal at high currents and 
the voltage drop is just ordinary Ohmic loss. There is a 
discontinuous change between the two regimes. The FF 
model is intended to simulate a situation in which the 
links are weak because flux pinning in them is weak. In 
this case, I c of the link represents the depinning current. 
The losses for I > I c are from the movement of flux along 
the boundary. Experimental measurements on individual 
grain boundaries show a cross-over between the two sorts 
of behaviour. Fig. || shows 7 — V characteristics for a 10 
degree grain boundary. At low applied field, the bound- 
ary resembles our JJ model, while at high field it is closer 
to the FF model. 

As to the more fundamental issue of whether our model 
is really valid for superconductors, we note that the main 
missing element is the phases of the grains. This may 
not be as bad an approximation as at first sight, how- 
ever. Our nonlinear network can be considered as an 
approximation to a true Josephson junction array in the 
overdamped limit. The J J model can represent a set of 
overdamped junctions with no pinning in which the Mc- 
Cumber parameter (3j = hR 2 /(2e/ c oC) 1 / 2 is less than 1. 
R and C are the resistance and capacitance of the junc- 
tions. Our output Itot then represents the DC component 
of the actual output of real overdamped junctions. The 
FF model is the case where pinning of Josephson vortices 
dominates the transport. 

The output of our calculations will be a Vtot(Itot) re- 
lation for the system as a whole. Thus, we feed a fixed 
current in at one end and collect it at the other, and mea- 
sure the voltage drop. We wish to compare this result to 
experimentally measured I — V characteristics such as 
that shown in Fig. |] and determine the extent to which 
the percolation is responsible for the macroscopic 
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FIG. 3. Experimentally measured I — V characteristics for 
a 10 degree YBCO grain boundary at various applied fields. 
At low fields, the I — V characteristic resembles our J J model. 
As the applied field increases, there is a crossover at 5.5 T to 
a behavior which is similar to our FF model. Data courtesy 
of N. Heinig and D.C. Larbalestier. 
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FIG. 4. Experimentally measured I — V characteristics for 
IBAD film at various applied fields. Again there is a crossover 
at roughly 5.5 T. As in our model (Fig. 9) the shape of the in- 
dividual grain boundaries I — V is preserved in the macrosopic 
system. Data courtesy of N. Heinig and D.C. Larbalestier. 

electrical properties and in particular, the critical cur- 
rent. In particular, we wish to understand whether the 
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presence of randomness and nonlinearity washes out most 
of the information about the individual resistors. 



III. ANALYTICAL APPROXIMATIONS 

In this section, we present an analytic approximation 
for the critical current of a binary lattice. Suppose we 
have a distribution of nonlinear resistors on an infinite 
lattice with each resistor possessing some critical current 
and occupation probability. (We shall take a lattice con- 
stant of unity to avoid distinguishing between currents 
and current densities.) This implies that the critical cur- 
rent can be found from the prescription: 

Ictot = mins'y]li t \ (1) 
ees 

where we take the minimum over all surfaces which sep- 
arate the electrodes p~2]. I are the wires which pierce 

S, and Ic is the critical current of wire I. This explicit 
but highly formal equation does not give any prescription 
for computation of I c tot except in very small systems. It 
says only that the network is only as strong as its weak- 
est hypersurface; in two dimensions, its weakest curve. 
However, because I ctot depends only on ijf , Eq. [l] tells 
us that the total critical current of the lattice depends 
only upon the critical currents of the individual elements 
and not upon their dissipative behaviour at / > I c . 

Consider two adjacent infinite surfaces S and S' , both 
of which separate the electrodes. The maximum super- 
current which can be carried through S, considered en- 
tirely by itself, is clearly 
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However, it is not clear that the current can penetrate 
from S to S 1 . Indeed we expect the critical current of 
the infinite lattice to be less than Ic°°\ which serves as 
an upper bound for I c t t, independent of model and di- 
mensionality. 

To obtain a better approximation for I c t t, consider a 
binary model, where there are only two types of resistors 
I c i and I C 2 with occupation probabilities p and q = 1 — p 
respectively. Since we are only interested in the critical 
current, the exact I — V characteristic of these resistors 
in unimportant. For this model the bound of Eqn. ^| is 
7* ot = pl c \ + qI C 2- As q increases, there will be a large 
increase in the the critical current when q > p c . This is 
obvious since for q > p c we get an infinite cluster of I C 2 
resistors running through the lattice, allowing much more 
current to be transported. However, it is not true that 
Ictot = Id for q < p c - This is most easily understood in 
the following manner. 

Consider an M x N square lattice and suppose there 
is a current / > NI C \ transported horizontally across the 
lattice. Due to the geometry of the lattice, we divide 



the resistors into two groups. We will call those resistors 
running perpendicular to the average current flow "row 
resistors" and term a line of these to be a "row". The 
rows are connected by what we will call "column resis- 
tors" and the resistors connecting two rows make up a 
"column", as shown in Fig. |^. The current flow can be 
considered a superposition of two current distributions, 
NI c i running directly through the lattice (along the col- 
umn resistors) , and a percolating current I — NI c i . It is 
this percolating current that we are interested in, so we 
can "subtract off" an I c \ resistor from each column resis- 
tor. A modified dilute resistor network lattice results and 
is shown in Fig. ^ The columns of this lattice have holes 
with probability p, and J C 3 = I C 2 — I& critical current 
resistors with probability q = 1 — p. This modified lat- 
tice has a critical current I c m The total critical current 
is then I ctot = I cl + I cM - 




FIG. 5. Picture of a two-dimensional network. The current 
runs to the right. We divide the resistors into those parallel 
to the current (column) and perpendicular (row). The actual 
lattice is then reduced to a dilute network by subtracting out 
a uniform current to the right. 



I c m is determined by two factors. First, the amount of 
current that can be transported across any given column 
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of the lattice, and second, whether or not that current can 
be redistributed to the J C 3 resistors in the next column. 
There are two different regimes of behavior distinguished 
by whether is less or greater than I c \ . The first case 
is simplest. 

When I C 3 < I c i (I c2 < 27 c i), it is fairly easy to re- 
distribute the current between columns. Any current 
coming across a column resistor can be shunted sideways 
along a row and redistributed to the next column as long 
as another column current doesn't get in the way. We can 
approximate the probability that a given column current 
can be redistributed along a row without being interfered 
with in the following fashion. It is equal to the probabil- 
ity that there is an J C 3 resistor directly across from the 
given column current, plus the sum of probabilities that 
there is a path at a site s steps along the row , and no 
paths either coming in or out of the row before that point. 
By only summing over steps to the right (or left) of the 
given column current, we avoid the problem of interfering 
with the paths of other incoming column currents. This 
redistribution probability can then be written as: 

oo 

p r =«x^(l- ? f = -^ (3) 

s=0 Lq q 

Since we have a volume fraction q of J C 3 row currents, 
we can then write down the averaged critical current den- 
sity as: 

htot = Id + {Ic2 - Id) x q xp r 

= hl + {Ic2-Id) X — j. (4) 

2q - q A 

It should be evident that this is only an approximate for- 
mula. The derivation neglects the possibility that the 
percolation current goes backward, for example. How- 
ever, we expect this path to make a smaller contribution 
than the ones calculated, since they must involve rela- 
tively rare configurations of resistors. 

A plot of Eq. [| for a 15 by 15 lattice is shown in Fig. 
^, where the results are compared to numerical data. It 
appears to overshoot the critical current by a substantial 
amount. We will see, however, that this is due to the 
finite size of the lattice used (15 by 15). For the infinite 
lattice case, the approximation should be very good if 

Ic3 < Id- 

If 2c3 > I c \ then this calculation will fail, since it is no 
longer as easy to shunt the current sideways along a row 
of resistors, and hence the redistribution probability will 
be different. A calculation similar to that above gives a 
different critical current: 

Ictot = Id + q{Id(Pr - Pr) + {Ic2 - Id)p' r } (5) 

where: 

oo 

p' r=(Z x^(A) s = -^ (6) 
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FIG. 6. Comparison of probabilistic theory to numeri- 
cal data. The critical current for a square lattice where 
Ici/Ic2 = 0.8 is approximated for varying occupation proba- 
bilities p and compared with the numerical data for a 15 x 15 
lattice. 



While this analyisis can clearly be extended to distri- 
butions of more than two resistors, it quickly becomes 
much more complicated. In addition, it tells us noth- 
ing about dissipative behaviour at currents greater than 
I c t t- Hence we now turn to numerics. 

IV. NUMERICAL CALCULATIONS 

Assume that we have a nonlinear resistor network and 
that the distribution of resistive elements is known. If a 
current is imposed across the network, a voltage drop re- 
sults. We wish to calculate this voltage drop for the entire 
lattice. The correct distribution of currents in the lat- 
tice can be found by solving Kirchoff's equations, which 
uniquely determine the current distribution in the net- 
work. In a network of linear resistors, this reduces to 
solving a set of linear equations, for which many stan- 
dard methods exist. For nonlinear resistors, we are not 
aware of any general method. The method of Leath and 
Tang |L0]], also used by Hinrichsen et at, pJ^ , works only 
for the JJ case. The approach we take is to search among 
all distributions of current which satisfy current conserva- 
tion and transport the imposed current across the lattice. 
We then determine which distribution in this class also 
has a unique voltage at each node. 

To search among all current distributions in this class, 
we begin by choosing a distribution which is thought to 
be close to the actual one. This choice must conserve cur- 
rent at each lattice node, and must also transport the im- 
posed current across the lattice. One possibility would be 
a distribution where the current moves uniformly across 
the lattice. This initial guess obviously needs to be modi- 
fied. We do this by superimposing circulation currents on 
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top of the initial distribution as shown in Fig. ^. While 
these circulation currents contribute nothing to the net 
transport of current, they alter the path that the cur- 
rent takes. Current conservation is obviously satisfied 
for all values of the circulation currents. The circulation 
currents have become our free variables. It is an elemen- 
tary result of homology theory JL4] that all circulating 
currents can be generated by the addition of currents 
circulating around plaquettes. Thus our search method 
can ultimately produce all possible currents. 



Circulation 
Current 




Initial 
Guess 

FIG. 7. Initial guess for the current distribution and impo- 
sition of circulation currents 



The current distribution satisfying the voltage Kirchoff 
law must still be found. We note that this requirement 
can be expressed in terms of a minimization principle. If, 
for a given resistor, we define the quantity: 



9j 



Vj(i)di 



(7) 



where Vj(i) is the voltage-current characteristic of the 
jth resistor (which we assume to be known) and Ij is 
the current flowing through the jth resistor, then solving 
Kirchoff 's laws is equivalent to minimizing: 



9j 



(8) 



where the sum is over all resistors in the lattice. We 
prove this by varying Q with respect to the currents such 
that current conservation is preserved. That is, vary Q 
with respect to a circulation current: 
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(9) 



Since the circulation current Ck travels around a specific 
loop, variations in Ck only affect the currents along that 
loop. In addition the effect is unit linear in Ck- 



SI 

— - — 1 (Jj- on loop k) 
dCk 



Hi. 

SC k 



(otherwise) 



where the sign is positive in the direction of Ck- 

sg_ 

SCk 



j£loop 



(10) 
Then, 

(11) 



This is valid for any loop in the lattice and hence we see 
that the minimization of Q with respect to the circulation 
currents is equivalent to solving Kirchoff 's voltage law. It 
is interesting to note that in the case of linear resistors, 
Q reduces to Q = (1/2) £\ Rjlj- So for the linear case, 
or indeed any network where all the voltages have the 
same power law dependence upon current, minimizing Q 
is equivalent to minimizing the power loss of the lattice, 
a familiar result. 

It seems unlikely to us that this variational principle 
for nonlinear systems is not known. However, we have 
not been able to find it in the literature. We note that 
this method is not at all limited to this electrical exam- 
ple. It should be applicable to nonlinear transport of any 
conserved quantity. Fluid flow in random porous media 
is one candidate. 

Due to its robustness and applicability to highly disor- 
dered systems with many variables we chose a simulated 
annealing algorithm to perform the minimization. Deter- 
mining whether an optimization algorithm has found the 
global minimum can be difficult. However, in the non- 
linear percolation case with monotonic I — V relations, a 
check can be made. At the global minimum, the voltage 
drop across the lattice should be a constant regardless of 
the path taken. By calculating the voltage drop along 
various paths, the accuracy of the solution can be deter- 
mined. 

The main advantage of this algorithm is its general- 
ity. It is immediately applicable to elements of arbitrary 
nonlinearity, as well as any geometrical arrangement of 
these elements. We also note that this algorithm sup- 
poses an imposed current across the lattice. An equiva- 
lent minimization principle for an applied voltage can be 
formulated by minimizing the sum: 



(12) 



where the sum is over the resistors, and the free variables 
are the voltages at each node Q. 

While we are mainly interested in the infinite lattice 
case, realistically we must work with lattices of finite size. 
While one generally approximates the infinite case by av- 
eraging over many small lattices, the finite size fluctua- 
tions can still be important. 



G 



To understand the effect finite size has on the critical 
current, take a binary model on a square lattice. Eq. |^ 
gave us the bound 



Id < 



Ir. 



N 



<plcl + (l-p)lc 



(13) 



for an infinite lattice. However, for a finite lattice, this 
upper bound decreases. For a N x N lattice of nodes, we 
have a resistor lattice with M = N — 1 columns and N 
resistors in each column Fig. g, We are interested in the 
fluctuations. The occupation probabilities may be p and 
q = 1 — p, but we are obviously going to see large fluc- 
tuations from the mean in a smaller lattice. So what we 
want is the expectation value for the maximum number 
of I c \ resistors in any of the M columns, since this will 
be the limiting factor on the critical current. In order to 
find this expectation value, we need the probability V(a) 
that the maximum number of I c ± resistors in any of the 
M = N — 1 columns is a (out of N) . 

To derive an expression for V(a), define /(a) to be the 
probability that a specific column has 'a' I c i's in it. Then 

/(«)= H ™ v P a (l- P ) N - a . (W) 
a! (TV — a)! 

Now define Pm{x,cl) as the probability that x columns 
(out of M) have a I c \ resistors in them. 



P N (x,a) 



M 



-/(af(l-/(a)) 



M-x 



(15) 



x\(M-x)V 
where < x < M. V(a) is then given by: 

Ha) = Ef=i [Pm (x, a) x FlL+i p m- x (0, i)\ a^N 



P(N) 



a = N 
(16) 



The above equation is simply the probability that x 
columns in M have a R\ resistors, multiplied by the prob- 
ability that there are no columns (in the M—x remaining) 
that have more than a R\ resistors in them, summed over 
x. 

We then use V(a) to find the expectation value of a: 



Finally we get the critical current for the lattice, 
al c i + (N- a)I c2 



Ir. 



N 



(17) 



(18) 



This analysis is only relevant for small lattices, since 
for very large lattices a will be very close to p . However, 
for smaller lattices this effect can indeed be the limiting 
factor, as can be seen by looking at Fig. || where we plot 
the critical current for a 15 by 15 lattice with I c i = 0.8 
and I C 2 = 1.0 for various occupation probabilities. The 
critical current is limited by the finite size of the lattice 
and that the critical current for the infinite case would 
be larger. 
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FIG. 8. Finite size calculation. Icxjlci = 0.8 The limiting 
effects of finite size on the critical current are calculated for 
a 15 x 15 square lattice and compared with the numerical 
results. 



V. I - V RELATIONS 

The current-voltage characteristics of both J J and FF 
binary models were studied. For the sake of computa- 
tional ease, the calculations were performed on square 
lattices, usually 10 by 10. The I — V characteristics from 
50 lattices were averaged to obtain the final result. It 
should be noted that the JJ I — V characteristic used 
in the calculations is slightly modified from that which 
was previously discussed. Instead of a strictly discontin- 
uous model, the transition at I c from V — to V = IR 
was made to be linear over a small current range 51. In 
addition to being more physically realistic, this modifi- 
cation made the convergence of the simulated annealing 
algorithm much quicker. 

The most important result of the numerical simula- 
tions is that the overall shape of the individual elements 
V-I characteristics is preserved when the elements are 
combined into a lattice. Fig. [)] compares a 10 by 10 lat- 
tice composed of JJ resistors with one of FF. p = 0.5, 
W-f C 2 = 0- 75 and R1/R2 = 0.75 in both cases. We 
see that the resulting overall characteristics are very dif- 
ferent. This in itself is not surprising as the individual 
components have different properties. However, the fact 
that the overall characteristics should resemble that of 
the individual elements so closely is surprising. Statis- 
tical averaging, even with strong nonlinearity, does not 
wash out the underlying input. 

Since both JJ and FF resistors are linear at high cur- 
rent, it is clear that the collective behavior of a lattice 
of resistors will be linear at high currents. In addition, 
the slope (conductivity) will be independent of the model 
used. We also expect that this high current behavior of 
the JJ and FF V-I curves is offset vertically, so that for 
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JJ we expect the voltage to go as V = ITZ and for FF, 
V = (X - ifjK. {if ^ I c tot necessarily ) This can be 
seen in Fig. |[ The nature of the transitions between the 
superconducting and the high current behavior as well as 
the current range over which the transition should occur 
is not as evident. 

2.5 | , , , , , , 1 
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FIG. 9. Comparison of the I — V characteristics for bi- 
nary JJ and FF models on a 10 by 10 percolative lattice. 
Both models possess the same distribution of microscopic crit- 
ical currents and normal state resistances. (I c i/Ic2 = 0.75 , 
R1/R2 = 0.75 ,p = 0.5) The shape of the individual circuit 
elements I — V characteristics is preserved when they are com- 
bined into a lattice. 

Our calculations show that the net transition behavior 
closely mimics that of the individual elements. Quali- 
taivcly, we believe the explanation is as follows. While 
Ictot is independent of the choice of model, as there is no 
difference in the behaviour of JJ and FF resistors below 
I c , the behaviour above I c tot is highly model dependent. 
FF resistors can conduct at all voltages, and hence the 
transition from superconductivity to dissipative conduc- 
tion for a FF lattice is smooth and gradual. However, 
the JJ resistors have an almost discontinuous jump in re- 
sistance at I c . It is harder to satisfy Kirchoff's laws at 
low voltages, since these low voltages are in effect "for- 
bidden" . Thus the rise in voltage at I cto t is very steep, 
and the breadth of the transition rather narrow. 

Exactly how narrow the transition to dissipative be- 
haviour is depends on the probability distribution of the 
critical currents. For the binary model, as I C 2~Id grows, 
the transition broadens. This is true for both JJ and FF, 
(although it is easiest to see in the J J case). We illustrate 
this in Figs. |lfj| and |TT] where the effects of varying I c i/I C 2 
for lattices of fixed occupation probability p(q = 1 — p) 
and fixed R1/R2 for both J J and FF lattices are shown. 
If on the other hand we vary the occupation probabil- 
ity p as in Figs. |l2| and [l3], we see that the transition 
from superconductivity to the high current behavior is 
broadest at p — 0.5, and gets much narrower towards 



p = or 1. Thus we see that as the probability dis- 
tribution for the critical currents of the individual ele- 
ments flattens out, the current range occupied by the 
transition from superconductor, to normal conductor be- 
comes much larger. This statement is still relevant for 
non-binary systems which have a distribution of I — V 
characteristics V(V(I)). 

Very detailed comparison between theory and experi- 
ment is beyond the scope of this paper. It would require 
the experimental determination of the entire individual 
resistor distribution as well as macroscopic I — V curves. 
However, we can get partial information about the former 
and compare to the latter as a preliminary exercise. 

Fig. U shows I — V characteristics measured across a 
single 10° grain boundary at various applied fields. Vary- 
ing the magnetic field applied to a grain boundary will 
change the grain boundary's properties. The critical cur- 
rent, in general, decreases with increased applied field. 
The high current resistance is dependent on the applied 
field as well. We however wish to focus for now on the 
shape of the boundary's I — V characteristic. At low 
fields the grain boundary displays a Josephson junction 
type behaviour which resembles our JJ model. At about 
6 T there is a transformation to a different high field be- 
haviour caused by the depinning of flux along the grain 
boundary. We modeled this high field behavior with our 
FF model. 

How does a sample consisting of many grain bound- 
aries behave? Fig. [j] gives the voltage current charac- 
teristics for an link comprised of many grain boundaries, 
made by ion-beam assisted deposition (IB AD). At low 
applied fields, the overall I — V characterisic looks some- 
what like a Josephson junction. As the field is increased 
the curvature of the characteristic changes at about 5 T, 
and it begins to look much more linear. It appears that 
the new dissipative behavior is dominated by flux creep 
as opposed to a Josephson type dissipation. Our resistor 
model supports this. We have shown through our resistor 
model that the microscopic single grain boundary (resis- 
tor) characteristics can (at least in the case of JJ and 
FF) be carried through by percolation to the macroscopic 
many grain boundary case. A macroscopic sample com- 
prised of grain boundaries acting as Josephson junctions 
will itself have a very sharply rising I—V characteristic, 
while a sample where the dissipation is dominated by flux 
creep will resemble its individual components as well. 

Comparison of our theory with the experiment sug- 
gests that the grain boundaries dominate the dissipative 
properties of the macroscopic superconductor, even in 
applied magnetic fields. If this is the case, then the ef- 
fects of the grains themselves are less important for the 
IBAD samples. 
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FIG. 10. 10 by 10 binary FF lattice. I - V characteristics 
for various I c i/I C 2 arc shown. p(I c i) = 0.5 , R1/R2 = 0.75 
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FIG. 11. 10 by 10 binary JJ lattice. I — V characteristic 
for various I c i/I C 2 are shown. p(I c i) = 0.5 , R1/R2 = 0.75 
There is substantial broadening of the transition near I c due 
to the finite lattice size. 



VI. CONCLUSION 

Conduction in polycrystalline high-T c superconductors 
is percolative, each percolative path being nonlinear. We 
believe that the model considered in this paper is general 
enough to capture the essence of this phenomenon. It re- 
quires powerful numerical techniques: the one introduced 
in this paper should be flexible enough to handle virtually 
any distribution of grain boundary conductances. For 
simple distributions, for example binary models, prob- 
abilistic reasoning can give insight into, and even good 
quantitative approximations for, the magnitude of the 
macroscopic critical current. 
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FIG. 12. 10 by 10 binary FF lattice for various occupation 
probabilities p(I c i)- I c i/I C 2 = 0.5 , R1/R2 = 0.75 
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FIG. 13. 10 by 10 binary JJ lattice for various occupation 
probabilities p(I c i)- Ici/Ic2 = 0.5 , R1/R2 = 0.75 Again, 
finite size transition broadening may be seen. 

Not only the critical current is of interest for applica- 
tions. One also wishes to know what the full I — V char- 
acteristic is telling us about the properties of individual 
grain boundaries. What emerges from the calculations 
is that the shape of macroscopic cur rent- volt age relation 
mirrors, to a remarkable extent, the shape of the under- 
lying current-voltage relations of the grain boundaries. 
This is counterintuitive because the latter have a distri- 
bution which may be broad, and one would expect this to 
produce a broadening which would wash out discontinu- 
ous behavior near I c tot- This does not happen. Nor does 
it happen that, for example, there is a two-step behavior 
for binary distributions. Instead, we find that the full 
I — V always has nearly the same shape as the individual 
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/ — V's, but with parameters which are averages over the 
underlying distribution. 

Our method allows us to find the full spatial distri- 
bution of current. We have not yet made use of this 
information. It is measurable by magneto-optical or Hall 
probe methods, and surely contains information about 
the spatial distribution of microscopic critical currents. 
This represents a possible future direction for this work. 
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